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The transition form factors of B( a i — * S, with S denoting a scalar meson, are investigated 
in the light-cone sum rules approach. The numerical values are approximately twice as large 
as that estimated in the light-front quark model and QCD sum rules approach. Using these 
form factors, we present the analysis of the decay rates for B —> ao(1450)i^(, B — » Kq (1430)//, 
B s — > J\o(1430)i^i and B s — ► /o(1500)ii with I — e,fj,,r. The results indicate that magnitudes of 
BR(Bo — > ao(1450)/i?i) and BR(B S — > Kq(1430)Wi) can arrive at the order of 10~ 4 , which can 
be measured in the future experiments to clarify the inner structure of scalar mesons. It is also 
observed that BR(B -> Kq (1430)r+r") and BR(B S / (1500)t + t~) are an order of magnitude 
smaller than the corresponding channels of e + e~ and fj> + fx~ final states due to the heavily suppressed 
phase space . Moreover, the longitudinal lepton polarization asymmetry for B — > ifo(1430)// and 
B s — » /o(1500)/r are also investigated, whose values are close to —1 for the e + e~ and fi + fj,~ pair 
except the region close to the end points. 

I. INTRODUCTION 

The inner structure of scalar mesons has been controversial for over three decades, which makes them one of the 
alluring issues in contemporary particle physics. In particular, the existence of some physical states, such as /o(1370), 
is still in dispute due to the absence of convincing evidence [1] . It is suggested that the scalar mesons with the masses 
below and above lGeV can be organized into two nonets in terms of their spectrum [2, 3]. The flavor singlet /o(600) 
(or cr), /o(980), the isodoublet A"q(800) (or k), and the isovector ao(980) constitute the nonet below 1 GeV; while 
/o(1370), /o(1500)// (1710), A' *(1430) and a (1450) form the other one near 1.5 GeV [2, 3]. Up to now, there is no 
general agreement on the nature of these states [4] due to the ambiguity existing in all the available interpretations 
including conventional qq states [3], glueball, hybrid states, molecule states [5-7] as well as tetra-quark states [8] and 
the superpositions of these contents [9-11]. Among all the scalar mesons, the A"q(1430) is predominantly viewed as 
su or sd state in almost all the models. Hence, it is justified to assume that ao(1450), Aq(1430) and /o(1370) being 
in the same nonent are respectively the ud, us and uu + dd states, based on the naive quark model, which is also the 
picture of scalar mesons adopted in this paper. 

Studies on the strong and electromagnetic decays of scalar decays have been received extensive interests in the 
literature [12-15]. Besides, the production properties of scalar mesons in ttN scattering, pp annihilation, 77 formation 
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and heavy meson decays can also serve as an ideal platform to explore the underlying structures of scalar mesons 
as well as the non-perturbative dynamics of QCD. Thanks to the progress of accelerator and detector techniques, 
both Belle and BaBar have observed strong indications of scalar mesons within a broad spectrum between 1.0 and 
1.5 GeV in B meson decays [16]. In this work, we will focus on the semi-leptonic weak production of scalars in 
-B( s ) decays, which are relatively clean compared to the hadronic decays from the theoretical viewpoint. Moreover, 
scmilcptonic decays of B^ mesons are also of great importance to determine the quark-flavor mixing matrix — the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [17, 18] and testing its unitarity under the requirement of the standard 
model (SM). 

The main job of investigating the semi-leptonic decays of B^ to the scalar mesons (S) is to properly evaluate 
the hadronic matrix elements for B^ — > S, namely the transition form factors, which arc governed by the non- 
perturbative QCD dynamics. Several methods exist in the literature to deal with this problem, such as simple quark 
model [19], light-front approach [20-22], QCD sum rules (QCDSR) [23, 24], light-cone QCD sum rules (LCSR) [25-27], 
perturbative QCD factorization approach [28-30]. The QCD sum rules approach, which is a fully relativistic approach 
and well rooted in quantum field theory, has made a tremendous success; however, short distance expansion fails in 
non-perturbative condensate when applying the three-point sum rules to the computations of form factors in the large 
momentum transfer or large mass limit of heavy meson decays. As a marriage of standard QCDSR technique [23, 24] 
and theory of hard exclusive process [31-38], LCSR cure the problem of QCDSR applying to the large momentum 
transfer by performing the operator product cxpansion(OPE) in terms of twist of the rcvclcnt operators rather than 
their dimension [39]. An important advantage of light-cone QCD sum rules is that it allows a systematic inclusion 
of both hard scattering effects and the soft contributions [40]. Phenomenologically, LCSR has been widely applied 
to investigate the semi-leptonic decays of heavy hadrons [41-44], radiative hadronic decays [45-47], non-leptonic two 
body decays of B meson [48-51] and strong coupling constants [52]. 

In the present work, wc would like to adopt LCSR approach to study the rare decay of B^ — > S. The essential 
inputs in the light-cone QCD sum rules is the hadronic distribution amplitudes other than vacuum condensates in 
the QCD sum rules. It is known that LCDAs are non-perturbative functions, which describes the hadronic structure 
in rare parton configurations with a fixed number of Fock components at small transverse separation in the infinite 
momentum frame. In an attempt to accommodate the experimental data, there have been continuous interests 
concentrating on the research of pre-asymptotic corrections to the distribution amplitudes of hadrons in the exclusive 
reactions over two decades. In particular, the leading twist and twist-3 distribution amplitudes of scalar mesons have 
been worked out in [3, 53] based on the QCD sum rules and conformal symmetry hidden in the QCD Lagrangian and 
we will use these amplitudes in this paper. 

The paper is organized as follows: In section II we present the effective Hamiltonian responsible for the b — > u, s 
transitions in the standard model, where the parameterizations of hadronic matrix elements are also collected here. 
Based on the trace formulae and equation of motion, we also derive the relations among form factors /+ (<j 2 ), f-(q 2 ) 
and /t(<? 2 ) in the large recoil and heavy quark limit. The Gegenauber moments of twist-2 and twist-3 distribution 
amplitudes obtained in the QCD sum rules are collected in section III. Then the sum rules for the various form 
factors on the light-cone are derived in section IV with the standard correlation function to the leading Fork state. 
After grouping the input parameters, the numerical computations of form factors in light-cone QCD sum rules are 
performed in section IV. Subsequently, we apply these form factors to analyze the decay rates of Bq — ► ao(1450)^z, 
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Bq — > Kq(1430)U, B s — > ^(1430)/^ and f? s — > /o(1500)Zi as well as the longitudinal lcpton polarization asymmetry 
for the modes induced by the flavor-changing neutral current, where a brief analysis on comparisons with the results 
obtained in the light-front quark model and QCD sum rules are also included in this section. The last section is 
devoted to the conclusions. 



II. EFFECTIVE HAMILTONIAN AND PARAMETERIZATIONS OF MATRIX ELEMENT 



A. Effective Hamiltonian for b — > u,s transition 



Integrating out the particles including top quark, W ± and Z bosons above scale fi = 0(rrib) , we arrive at the 
effective Hamiltonian responsible for the b — > u transition 

G 

H eff (b -> ulvi) = -^V ub u ltl (l - 75 )^V(1 ~ 75>l + h.c. , (1) 

where V v b is the corresponding Cabbibo-Kobayashi-Maskawa (CKM) matrix clement and / = (e, /i, r) . 

Similarly, the effective Hamiltonian revelent to the flavor-changing neutral current (FCNC) transition b — > s can be 
derived as 

H eff (b - sll) - ^ Vtb V t :[Cf f ( M ) s 7 m(1 " W - Js)l + C 10 s^(l - l5 )bl^ l5 )l 

_ 2m b C 7 {») lh)(fh - m] + h. c . , (2) 

T 

where we have neglected the terms proportional to V u bV* a on account of IKibKTs/^tfc^tsI < 0-02 . The C, involved 
in Eq. (2) are the Wilson coefficients and their particular expressions are given in Ref. [54]. We should emphasize 
that the Wilson coefficient Cio does not renormalize under QCD corrections and hence is independent on the energy 
scale fi ~ 0(m b ), since the operator Oiq = S7 M (1 — 75)^7^75^ can not be induced by the insertion of four-quark 
operators due to the absence of Z boson in the effective theory. Moreover, the above quark decay amplitude can also 
receive additional contributions from the matrix element of four-quark operators, which are usually absorbed into the 
effective Wilson coefficient Cg (/i). To be more specific, we can decompose Cg (n) into the following three parts 
[55-61] 

C 9 e// (/i) = Caift) + Y SD (z, s') + Y LD (z, s'), (3) 

where the parameters z and s' are defined as z = m c jm bl s' = q 2 /m b . Ysd{z,s') describes the short-distance 
contributions from four-quark operators far away form the cc resonance regions, which can be calculated reliably 
in perturbative theory. The long-distance contributions Yld{z,s') from four-quark operators near the cc resonance 
cannot be calculated from first principles of QCD and are usually parameterized in the form of a phenomcnological 
Brcit-Wigncr formula, which will be neglected in this work due to the absence of experimental data on B( s ) — > J/ipS. 
The manifest expressions for Ysd{z, s') can be written as [54] 

Y SD {z,s') = h{z, 80(3(7! (m) + C 2 (m) + 3C 3 (m)+C4(m)+3C 5 (m) + C 6 (m)) 
~h(l, 8'){4C 3 (ji) + iC 4 (fi) + 3C 5 (/i) + C 8 0*)) 

-^(0, s')(C a (fi) + 3C 4 (/x)) + ^(3C 3 ( M ) + C4p) + 3C 6 (/i) + C 6 ( M )), (4) 
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with 



h(z,s') = 
h(0,s') = 



-Inz 



± + 1,- §(> + ,),! -^v 



Ill 



2 arctan ■ 



in for x = 4z 2 / V < 1 
for x = 4z 2 /s' > 1 



— —In — - — —Ins' 



(5) 



27 9 n 9 9 

Besides, the non-factorizablc effects from the charm quark loop can bring about further corrections to the radiative 
b — ► S7 transition, which can be absorbed into the effective Wilson coefficient as usual [62-65] . Specifically, the 

Wilson coefficient is given by [66] 

c 7 e// ( M ) = c 7 (^) + cu», 



with 



Gi(x) 



ia 8 [-?7 14/23 (Gi(^) - 0.1687) - 0.03C 2 (/x)], 
y 

x(x 2 -5x~2) 3x 2 \n 2 x 



(x - l) a 



4(x- l) 4 ' 



(6) 

(7) 
(8) 



where r\ = a s (mff)/a s (/i), x± = ra 2 /m 2 V) C^_ >S7 is the absorptive part for the b — > sec — > ,57 rescattering and we 
have dropped out the tiny contributions proportion to CKM sector V u bV* B . 

B. Parameterizations of hadronic matrix element 

With the free quark decay amplitude available, we can proceed to calculate the decay amplitudes for semi-leptonic 
decays of B q / — > S at hadronic level, which can be obtained by sandwiching the free quark amplitudes between the 
initial and final meson states. Consequently, the following two hadronic matrix elements 

(S(p)\s W b\B q ,(p + q)), (S(p)\8a ia ,TB<fb\B g >(p + q)) (9) 

need to be computed as can be observed from Eqs. (1) and (2). The contributions from vector and tensor types of 
transitions vanish due to parity conservations which is the property of strong interactions. Generally, the above two 
matrix elements can be parameterized in terms of a series of form factors as 



(S(p)\sw 5 b\B q ,(p + q)) 
{S(p)\sa llu j 5 q v b\B g ,(p + q)) 



-i\f+(tf)PM + f-tf)q»], 

(2p + q) M q 2 - (m| - m%) q^ f T {q 2 ) 



(10) 
(11) 



tub + ms 

Utilizing the covariant trace formalism introduced in [67] , the form factors at large recoil should satisfy the following 
relations 

2m B 



Mi 2 ) = 



-Mq% 



m B + m s 

where the corrections due to hard gluon exchange are neglected 
and f T (q 2 ) as 



/-(9 2 ) = 0, (12) 
|. We can also derive the relation 1 between f-(q 2 ) 



f T (q 2 ) = - 



m b - m q2 
ms - ms ' 



Mi 2 ), 



(13) 



1 This relation has been derived in Ref. [69], where the convention cr M „ = — 5 ("//Ok — 7^7n) is adopted . Therefore, the relations between 
f—(q 2 ) and /t(<? 2 ) there differ from that given in this paper with a minus sign. 
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with the help of the equation of motion in the heavy quark limit [70] . 



III. DISTRIBUTION AMPLITUDES OF SCALAR MESONS 



(S{p)\q 2 (x) lll q 1 (y)\0) = p u / fa^ I+i ")$ s ( V ), 

Jo 

(S(p)\q 2 (x) qi (y)\0) = m s [ due^* + ^<5> s (u,fi). 



In this section, we would like to collect some revelent information on the distribution amplitudes for scalar mesons, 
which are the essential ingredients in the sum rules for the B — > S transition form factors. 

Up to the leading Fock states, the light-cone distributions of scalar mesons made up of (72Q1 can be defined as [3, 53]: 



(14) 



where z = x — y, ms is the mass of corresponding scalar meson, u = 1 — u and u is the momentum fraction carried 
by the quark q 2 in the scalar meson. Here <£>g (u, fi) is of twist-2, <&|. (u, fi) and $g (it, fi) are of twist-3. <E>s (u, fi) and 
($5 (u, (j,), <I>g (u, fi)) are anti-symmetric and symmetric under the replacement of u — > 1 — u in the SU(3) limit owing 
to the conservation of G parity. To be more specific, the normalizations of $5 (it, fi), & s (u, fx) and (u, fi) are given 
by 



<S (p) |«j (ar) <r^gi 0) = -m s ( Pfl z v - PuZfl ) due^ x+ ^<S> s (u, fi) 



du$ s (u, fi) = fs, [ du<f> s s (u, fi)= [ du$ s (u, ft) = /<?. 
Jo Jo 

The vector current decay constant fs defined by 

(5(p)|g 2 7 M 3i|0) = /s^ 

should vanish in the SU(3) limit and can be related to the scalar density decay constant fs determined by 

(S (p) \q~2qi\ 0) = m s fs 

with the help of the equations of motion 

fs = fJ'sfs, fJ-s 



m s 



(15) 



(16) 



(17) 



(18) 



m 2 (fi) - mi (fi) ' 

with mi and m 2 being the masses of quarks q\ and q 2 respectively. It should be emphasized that scalar density meson 
decay constant fs depends on the renormalization scale fi, whereas the the vector current decay constant fs does not 
renormalize under the QCD corrections due to the conservation of vector current. 

In view of the conformal symmetry hidden in the QCD Lagrangian, the distribution amplitudes of scalar mesons 
"3?S (w, fi), $5 (u, fi) and <5> s (u, fi) can be expanded in the Hilbert space by Jacobbi polynomials with increasing 
conformal spin as 



®s(u,fi) = f s (fi)6uu 



B Q (fi) + J2 B m (fi) C 3 J 2 (2u - 1) 



= fs(jt) 



l+5> ro (/i)CV 2 ( 2tt -l) 



m= 1 



(19) 



$g (u,fx) = fs(fi)Guu 



1+ £& m (M)C r 3 / 2 (2«-l) 



rn—1 
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TABLE I: Decay constants and Gegenbauer moments for the twist-2 distribution amplitude $s of scalar mesons at the scale 
H = lGeV [3]. 



state /(MeV) 


Si 


B 3 


a„(1450) 460 ±50 


-0.58 ±0.12 


-0.49 ±0.15 


A" *(1430) 445 ±50 


-0.57 ±0.13 


-0.42 ±0.22 


/o(1500) 490 ±50 


-0.48 ±0.11 


-0.37 ±0.20 



TABLE II: Gegenbauer moments for the twist-3 distribution amplitudes $| and $5 of scalar mesons at the scale n — lGeV 
[53]. 



state 


ai(xl0~ 2 ) 


(1-2 




04 


bi(xl0~ 2 ) 




b 2 






ao 





-0.33 ~ 


-0.18 


-0.11 ~ 0.39 





<- 


< 0.058 


0.070 ~ 


0.20 


Ko 


1.8 ~ 4.2 


-0.33 ~ 


-0.025 




3.7 ~ 5.5 


r 


- 0.15 






fo 





-0.33 ~ 


-0.18 


0.28 ~ 0.79 





-0.15 


~ -0.088 


0.044 ~ 


0.16 



3/2 1x3/2 

where Gegenbauer polynomial C m (x) can be considered as a special type of Jacobbi polynomials P„l + 2( x ) ~ ^™ ( x )- 
Combining Eqs.(15), (18) and (19), the zeroth Gegenbauer moment B (fJ.) for twist-2 distribution amplitude $s(u,//) 
is given by 

B = Ms 1 - (20) 

Moreover, decay constants of scalar mesons and various Gengauber moments B m , a m and b m for both twist-2 and 
twist-3 LCDAs have been computed in Refs. [3, 53] based on QCD sum rules approach, which are collected in Table 
I and II. 

IV. LIGHT CONE SUM RULES FOR FORM FACTORS 

With the LCDAs of scalar mesons available, we are now in a position to derive the sum rules of transition form 
factors which are responsible for — > S decays. The basics object in LCSR approach is the correlation function 
in which one of the hadron is represented by the interpolating current with proper quantum number, such as spin, 
isospin, (charge) parity and so on; and the other is described by its vector state manifestly. Information on the hadronic 
transition form factor can be extracted by matching the Green function calculated in two different representations, i.e., 
phcnomenological and theoretical forms, with the help of dispersion relation under the assumption of quark-hadron 
duality. 
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A. Light-cone sum rules for the form factors f+(q 2 ) and f-{q 2 ) 

Following the standard procedure of sum rules, we consider the correlation function associating with the form 
factors / + (g 2 ) and /_ (g 2 ) determined by the matrix element 

n^(p,?) = - / d A xe^ {S {p)\T {j 2li {x), ji(0)}|0), (21) 



where the current j2^{x) = q 2 (x) 7^756 (x) describes the weak transition of b to 52 and ji(0) = b (0) 175(71 (0) represents 
the B qi channel. In addition, the vacuum-to-meson matrix element for the interpolating current can be given by 

2 

(B qi \bi l5 q\0) = T^—f Bqi . (22) 

1 1 m b + m qi 51 

Inserting the complete set of states between the currents in Eq. (21) with the same quantum numbers as B qi , we can 
arrive at the hadronic representation of the correlator (21): 

. (S (p) \q 2 (0) W b (0)| B qi (p + q))(B qi (p + q)\b (0) i7 59i (0) |0) 



n M (p, q) 



2 



m B n - + l) 

V- ■ (S (p) \q 2 (0) 7/,756 (0)| hjjp + q)) (h(p + q)\b (0) ?75gi (0) |0) 
^ m 2 -(p + g) 2 



where we have separated the contributions from the ground state and higher states corresponding to the B qi meson 
channel. Combining the Eqs. (10) , (22) and (23), the phenomcnological representations of correlation function (21) 
can be derived as 

TT ( \ m2 B q jB n 2 2 p ^(S^K+^M 2 )^ 

n^fe? = 7 -r^ ; 717 /+ (9 K + /-0? <?m + B ds 1 t ^ ' 24 ) 

m b + m 9l mi - p + g 2 ^0 w « - (P + 9) 2 



where we have expressed the contributions from higher states of the B qi channel in the form of dispersion integral 
with Sq 91 being the threshold parameter corresponding to the B qi channel. 

On the theoretical side, the correlation function (21) can also be computed in the perturbativc theory with the help 
of OPE technique at the deep Euclidean region p 2 , g 2 = -Q 2 < 0: 

IL>,g) = U^ CD (q 2 ,(p + q) 2 )p^ + n QCD (q 2 ,(p + q) 2 ) q ^ 

l lmllf^W) p Ihnljj^W) 
ds ■; — ■ — Pf, + / ds -; — ■ — 7^— g M . (25) 



Making use of the quark-hadron duality assumption 



P 1(s,q 2 ) = -lmnf CD ( S ,q 2 )Q(s-s h ), (26) 

7T 



with i = "+, — " and performing the Borel transformation 



-(P+9) 2 ,>1 

-(p+ 9 )V"=A/ 2 



with variable (p + q) 2 to both two representations of the correlation function, we can finally derive the sum rules for 
the form factors 



- / lmn^(.s, 9 2 )exp( -^f— )ds. (28) 
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b 

FIG. 1: The tree level contributions to the correlation function Eq. (21), where the current ji(0) represents the B qi channel 
and the current ,72(2;) describes the b — > 52 transition. 



To the leading order of a s , the correlation function can be calculated by contracting the bottom quark field in Eq 
(21) and inserting the free b quark propagator 

d 4 k e t(q -^ x 



(p, q) 



(2tt) 4 m 2 b - k 2 



(S (p) \q 2 (x) 7^75 (# + m b ) i-f 5 qi (0)| 0) , 



(29) 



which can be represented by Fig. (1) intuitively. It should be pointed out that the full quark propagator also receives 
corrections from the background field [71, 72] and can be written as 



(0|T{6 I (x)6,(0)}|0) = Sn 



d 4 fc 
(2^ 



}t — m& 
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d*k 



-ika. 



(30) 



where the first term is the free-quark propagator and = G^T°- with Tr[T a T b ] = \S ab . Substituting the second 
term proportional to the gluon field strength into the correlation function can result in the distribution amplitudes 
corresponding to the higher Fock states of scalar mesons. It is expected that such corrections associating with the 
LCD As of higher Fock states do not play any significant role in the sum rules for transition form factors [73] , and so 
can be safely neglected. 

Substituting Eq. (14) into Eq. (29) and performing the the integral in the coordinate space, we can achieve the 
correlation function in the momentum representation at the quark level as 



n M (p, q) = P M / du 



m l ~ (<? + u pY 



du 



-m b $ s (u) + itm 5 $| (u) + im 5 $5 (u) 



ml — u 2 p 2 + q 2 
m b - (<? + u vf 



m l - (q + u p) 2 



QCD, 2 



m b + u P ~ 1 

m l - {q + u p f 



(q z ,(p + qr)p» + n^ (q,(p + qr)q» 



(31) 



Combining the Eqs. (28) and (31), we can finally arrive at the sum rules for form factors f+(q 2 ) and f-(q 2 ) as 
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below 



f (<z 2 ) = 



with 



(m b + m qi ) 

tS qi J a q\ 



■ exp 



( m% 



du 



«0 

1 rn s 
uAP 6 



■ cxp 



$5 («) (r 



m 2 + imp 2 — ug 2 



2„2 



m b $ s (u) + m s ( u<P s s (u) + -$5 (it) 



(m b + m qi ) 



exp 



m B 

M 2 



du 



■ exp 



u 



u 

_L rns_ 

l 2 M 2 6 



-(so 



u p 



m 2 + uup 2 — uq 2 



Mcxp 



uM 2 

$5 («) (to 2 + u 2 p 2 - q 2 ) 



' M 2 J 
1 



1 - u p 2 + q 2 

2 + u 2 p 2 - q 2 



ms ( $| (u) + — $5 («) 



TOg 



6u 



K)exp(-^) 



(32) 



(33) 



P 2 ) + V( s o ~ I 2 ~ P 2 ) 2 + 4p 2 (mjj - q 2 ) 



2p 2 



(34) 



As can be observed from the sum rules (32) and (33), both twist-2 and twist-2 distribution amplitudes of scalar 
mesons can contribute to the form factor /+(g 2 ), whereas the other one /_(g 2 ) can only receive the contributions 
form twist-3 LCDAs and should be heavily suppressed in the large recoil region, which is also in agreement with the 
relations (12) presented in [68]. 



B. Light-cone sum rules for the form factor fr(q 2 ) 

As for the form factor fx (<? 2 ) involved in b — ► s transition, we start with the following correlation function 

fi M (p,?) = -J^xe^iS^lT^x),^)}^) (35) 
where the current J2^{x) is given by 

32n(x) = q2{x)a fll/ a"j 5 b(x) . (36) 

One can write the phenomenological representation of the correlation function at the hadronic level simply by 
repeating the procedure given above as 

Tor fg 

n„(p,?) = 2 q \ - . 2u -[(2p + q),q 2 -q,(m 2 B - m 2 )}f T (q 2 ) 

(m b + m qi ) [m% n - (p + q) \(m B + m s j 

ds 1 , 2 i-P^ 2 + QuiP ■ q)]Pt(s, Q 2 )- (37) 

i s - (p + q) 2 

On the other hand, the correlation function at the quark level can be calculated in the framework of perturbative 
theory to the leading order of a s as 

ft, (p, q) = [-JM 2 + ?m(p • q)\ f du^-^ 3 Us («) - mb ™ S 2 niu) A ■ (38) 

Jo ml - (q + up) [ o m l — (q + up) J 

Matching the correlation function obtained in the two different representations and performing the Borcl transfor- 
mation with respect to the variable (p + q) 2 , we can achieve the sum rules for the form factor /t(<Z 2 ) 



, 2 \ {m b + m qi ) (m B + m s ) ( m 2 B \ ( 1 f 1 du 
h GO = J-^I ex P ( jfi ) \ - 2 / ~ 6XP 



no U 



to 2 + uup 2 — uq 2 
vT 



to b tos 

$c (it) (u) 

6 V ' 3uM 2 5 V ' 



-^^ s (no) exp (*) L j I. (39) 

6 S ^ ; V Af 2 / m 2 + MgP — q 
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V. NUMERICAL ANALYSIS OF TRANSITION FORM FACTORS 

Now we are going to analyze the sum rules for the form factors numerically. Firstly, we collect the input parameters 
used in this paper as below [74-81]: 

Gp = 1.166 x 10- 2 GcV~ 2 , \V ub \ = 3.96±g;g| x 10~ 3 , 

\V tb \ = 0.9991, \V U \ = 41.61±g;|g x lO" 3 , 

m b = (4.68 ± 0.03)GeV, m s (lGeV) = 142MeV, 

m„(lGeV) = 2.8MeV, m d (lGeV) = 6.8MeV, 

mg = 5.279GeV, m Bs = 5.368GcV, 

f Bo = (0.19 ± 0.02)GcV, f Bs = (0.23 ± 0.02)GeV. 

It is noted that the input values for f B and f Bs are in agreement with the unquenched lattice results [79] fg = 
0.216 ± 0.022 GeV and f B , = 0.259 ± 0.032 GcV, and with the results from the QCD sum rules [80, 81]. The 
threshold parameter s can be determined by the condition that the sum rules should take on the best stability 
in the allowed Borel region. Besides, the values of threshold parameter should be around the mass square of the 
corresponding first excited state, hence they are also chosen the same as that in the usual two-point QCD sum rules. 
The standard value of the threshold in the X channel is so x = ( m x + Ax) 2 , where Ax is approximately taken to 
be 0.6GeV in the literature [82-85] . To be more specific, we adopt the threshold parameters Sq° = (35 ± 2)GeV 2 
and Sq" = (36 ± 2)GeV 2 corresponding to Bq and B s channels respectively, for the error estimate in the numerical 
analysis. 

With all the parameters, we can proceed to compute the numerical values of the form factors. In principle, the form 
factors f+(q 2 ), f-(q 2 ) and /t(<? 2 ) should not depend on the Borel mass M 2 in a complete theory. However, as we 
truncate the operator product expansion up to the leading conformal spin of distribution amplitudes for scalar mesons 
in the leading Fock configuration and keep the perturbative expansion in a s to leading order, a manifest dependence 
of the form factors on the Borel parameter M 2 would emerge in practice. Therefore, one should look for a working 
"window" , where the results only mildly vary with respect to the Borel mass, to make the truncation reasonable and 
acceptable. 

Firstly, we concentrate on the form factors at zero momentum transfer. As for the form factors /+(0) involved 
in Bq — > ao(1450)/^, we require that the contributions from the higher excited resonances and continuum states 
hold the fraction less than 20 % in the total sum rules and the value of /+(0) does not vary drastically within the 
selected region for the Borel mass. In view of these considerations, the Borel parameter M 2 should not be too large in 
order to ensure that the contributions from the higher states are exponentially damped as can be observed from Eqs. 
(32), (33) and (39) and the global quark-hadron duality is satisfactory; on the other hand, the number of Borel mass 
also could not be too small for the sake of validity of OPE near the light-cone for the correlation function in deep 
Euclidean region, since the contributions of higher twist distribution amplitudes amount the higher order of 1/M 2 to 
the perturbative part. Subsequently, we indeed find the Borel platform M 2 £ [10, 15] GeV 2 with the selected threshold 
parameter Sq° = 35GeV 2 as shown in Fig. 2, where the ratio of twist-3 distribution amplitudes in the total sum rules 
are also included for a comparison. Following the same methods, we can also further evaluate all the form factors 
/+(0), /_(0) and / T (0) associating with the decay modes B -> a (1450)lu h B -> K$ (1430)11, B s -> K$(1430)lvi, 
and B s — > /o(1500)ZZ, whose results have been collected in Table III- VI, where we have combined the uncertainties 
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FIG. 2: The solid line denotes dependence of form factors /+ at q 2 = responsible for the decay of Bo — > ao(1450)Z^ on the 
Borel window M% £ [10.0, 15.0] GeV 2 with the chosen threshold parameter s^ qi = 35GeV 2 . The dashed and the dot-dashed 
lines are the ratio of contribution from twist-3 distribution amplitudes and higher states of B qi channel in the total sum rules 
respectively. 



from the variation of Borel parameters, fluctuation of threshold value, errors of b quark mass, corrections from decay 
constants of the involved mesons as well as uncertainties from the Gengenbauer moments in the distribution amplitudes 
of scalar mesons. It can be observed that the errors on the form factors are estimated within the level of 20 % as 
expected by the general understanding of the theoretical framework. 

Next, we can further investigate the q 2 dependence of the form factors f + (q 2 ), f-{q 2 ) and /t(<? 2 ) based on the sum 
rules given in Eqs. (32), (33) and (39). One usually parameterize the form factors fi(q 2 )(i = +, —,T) in cither the 
single pole form 

Mi 2 ) = i — wt-> ( 41 ) 

1 - a l q 2 jm l B 



or the double-pole form 



/i(0) 

1 - a iq 2 /m 2 B + hq^/mi 



Ml 2 ) = , T^JZP^TJZT, , (42) 



in the whole kinematical region < q 2 < {in-B qi — nis) 2 , while non-perturbative parameters a% and hi can be fixed 
by the magnitudes of form factors corresponding to the small and intermediate momentum transfer calculated in the 
LCSR approach. Our results for the parameters at, hi accounting for the q 2 dependence of form factors /+, /_ and 
fx are grouped in Table III- VI, where the values estimated in other works are also given for a comparison. 
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TABLE III: Numerical results for the parameters fi{0), a,i and bi involved in the (single) double-pole fit of form factors (41), 
(42) responsible for Bo — > ao(1450)Z^ decay up to the twist-3 distribution amplitudes of scalar mesons, where the numbers 
derived in the covariant light-front quark model [86] are also collected for a comparison. 





h(0) 


a% 


h 


/+ 


1 04+ 020 

1 - u ^-0.20 


0.981°;°° 






0.52 [86] 


1.57 [86] 


0.70 [86] 


/- 


0.077«;^ 4 


1 ^9 +0 ' 07 




h 


0.66t >;- 


0.88l°,; ° 





TABLE IV: Numerical results for the parameters fi(0), a,i and bi involved in the (single) double-pole fit of form factors (41), 
(42) responsible for Bo — * Kg (1430)11 decay up to the twist-3 distribution amplitudes of scalar mesons, where the numbers 
derived in the covariant light-front quark model [87] and QCD sum rules approach [69] are also collected for a comparison. 





HQ) 


<h 




/+ 


07+0-20 
u - 3 ' -0.20 


0.86l° ;i 9 8 






0.52 [87] 


1.36 [87] 


0.86 [87] 




0.62 ± 0.16 [69] 


0.81 [69] 


-0.21 [69] 


f- 


0.0731°,;°,* 


2.50l°;f 7 


i o 9 +0.69 


h 


o.eoiS:iS 


0.69+°;?,? 






0.34 [87] 


1.64 [87] 


1.72 [87] 




0.26 ± 0.07 [69] 


0.41 [69] 


-0.32 [69] 



VI. DECAY RATE AND POLARIZATION ASYMMETRY 



With the transition form factors derived, one can proceed to perform the calculations on some interesting observables 
in phenomenology, such as decay rate, polarization asymmetry. In particular, the forward-backward asymmetry for 
the decay modes B — > A"q(1430)(7 and B s — > / (1500)2Z is exactly equal to zero in the SM [89, 90] due to the absence 
of scalar-type coupling between the lepton pair, which serve as a valuable ground to test the SM precisely as well as 
bound its extensions stringently. 

The semi-leptonic decay B — > ifg (1430)i7 is induced by flavor-changing neutral current. The differential decay 

TABLE V: Numerical results for the parameters /i(0), cn and bi involved in the (single) double-pole fit of form factors (41), 
(42) responsible for B a — ► ifo(1430)/^( decay up to the twist-3 distribution amplitudes of scalar mesons, where the numbers 
derived in the QCD sum rules approach [88] are also collected for a comparison. 





Mo) 


n, 


bi 


/+ 


0.83+°,;?* 
0.48 ± 0.20 [88] 


0.93i° :3 
1.25lg;S [88] 




/- 


0.07l±g;g 


2 4fi+ - 36 


-1 7 9+0.59 
^- lz -0.6i 


h 


S2+ 018 

u - dz -0.08 


77+ - 10 

u - ' ' -0.07 
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TABLE VI: Numerical results for the parameters fi(0), a,i and bi involved in the (single) double-pole fit of form factors (41), 
(42) responsible for B s — > /o(1500)££ decay up to the twist-3 distribution amplitudes of scalar mesons. 





Ho) 


a. 


bi 


f+ 


0M±° f 5 


-, 17 +0.06 
L - L '-0.05 




f- 


0.056t°;™ 


1 94+ 048 


0.52±?; 8 6 9 


h 


0.56l° ;iS 


i nq+ 08 





width of Bq — ► /Vq(1430)^ in the rest frame of -Bo meson can be written as [74] 

dT{B Q -> A' * (1430)2/) 1 1 



(43) 



dq 2 ( 27r )3 32mSo 

where u = (pk*(i430) + Pi) 2 and q 2 = (pi +Pi) 2 ', Px*(i430)7 Pi and Pi are the four-momenta vectors of Aq(1430), I and 
I respectively; |-^g ^ifj(i430)iz~| 2 i s the squared decay amplitude after integrating over the angle between the I and 
A"q(1430) baryon. The upper and lower limits of u are given by 



(E 
(E 



if *(1430) 



ff*(1430) ""if*(1430) 



- jEf - m 2 f 



K *(1430) 



Ety-(JE* K i 



J /f*(1430) '""if* (1430) 



(44) 



where E 



K *(1430) 



and Ei are the energies of Aq (1430) and I in the rest frame of lepton pair and can be determined as 



E 



1 ? 9 

m% — m, 

Ba 



AT *(1430) 



1 K* (1430) 1 



Ef 



1/2 



rfT (B q , -+ Sll) = G 2 F \V tb V ts \ 2 m%a 2 em A _ 4rA ■ 1/2 
ds' 1536tt 5 I s' / ^ s 



2r; . 

H T j a s + rid s 



where 



(45) 



2yV Wl 2 
Collecting everything together, we can arrive at the general expression of differential decay rate for B q i — > Sll as [87 



(46) 



2 i'i 1 / 1 1 1 1 

q /m B , n=mi/m B , r s = m s /m B , 



(l-r s ) 2 -2s{l + r s )+s 2 , 



a s = Lps 



a 



; ///+(g 2 ) n C 7 f T (q 2 ) 



l + \A\s 



/+ (v 2 ) 



6\C 10 f { [2(1 + r s )-s] 



/+ (? 2 ) 



(l-r s ) Re 



hid 2 ) (f-(q 2 ) 



U (<z 2 ) 



f-(q 2 ) 



The invariant dilepton mass distribution for Bq —> Aq(1430)Z2 as the functions of squared momentum transfer q 2 
are presented in Fig. 3. In the same way, we can also estimate the decay rates of Bq — > Aq(1430)ZL B s — > A'q (1430)ZPj 
and B s — > /o(1500)ZZ with I = e,fx,T based on the form factors calculated in light-cone sum rules. The results of the 
total decay width corresponding to these decay modes are grouped in Table VII, where the results obtained in other 
frameworks are also presented for a comparison. As can be observed, the decay rates of the electron- and muon- pair 
final states are practically the same, while the decay rate of tauon-pair channel is much smaller due to the heavily 
suppressed phase space. 
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FIG. 3: The invariant dilepton mass distributions for B a -> K£(U30)e + e~ , B -> A^(1430W and B -> A^(1430)t + t~ 
as functions of squared momentum transfer q 2 based on light-cone QCD sum rules. 



Another interesting observable in the decay of Bq — ► Kq (1430)11 is the polarization asymmetry of the final state 
charged leptons, which is very helpful to extract the information on the spin of them. The four-spin vector s M of a 
lepton can be defined in its rest frame as 



(S^r.s. = (0, C). 



(47) 
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TABLE VII: Numerical results for the total decay width of B -> a (1450)Zi/j, B -> #£(1430)//, B s -> K(j[(1430)ii/j and 
B s — * /o (1500)// with I = e, /i, r in the light-cone sum rules approach, together with the numbers estimated in QCD sum rules 
[69, 88] and light-front quark model [87]. 

Bo -> ao(1450)ei/e Bo K$ (1430)e + e~ S s #5(1430)ePe B 3 /o(1500)e+e _ 

LCSR l-8±S:i x 10 ~ 4 5.7t|* x 1(T 7 l-3io.I x KT 4 5.3±?j x 10~ 7 

LFQM 1.63 x 10~ 7 [87] 

QCDSR (2.09 ~ 2.68) x 10~ 7 [69] 3.6±|;f x 10~ 5 [88] 

Bo -> o (1450)/n^ Bo -» A" *(1430W B s -> KJ(1430)^ B s -> / (1500)/x+m" 

LCSR 1.8±g;? x 10~ 4 5.6±1 ;^ x 10~ 7 1-3±J:2 x KT 4 5.2^ j x 10" 7 

LFQM 1.62 x 10~ 7 [87] 

QCDSR (2.07 ~ 2.66) x 10~ 7 [69] 

Bo -> a (1450)rJ/ T B -> K$ (1430)t+t" B b -> Kq (1430)tP t B s -> / (1500)r+r" 

LCSR 6.3t| 4 x 10~ 5 9-8ls 2 5 4 x 10" 9 5.2±f;^ x 10" 5 1.2±gJ x 10~ 8 

LFQM 2.86 x 10" 9 [87] 

QCDSR (1.70 ~ 2.20) x 10~ 9 [69] 



The unit vector along the longitudinal direction of the lepton polarization is given by 



Pi 

\pi\ 



In this work, we mainly concentrate on the longitudinal lepton polarization asymmetry that can be defined as 

dr(e L £=i) _ dr(e L £=-i) 
Pl{ S ') = ' 



dr(e L £=i) dr(e L £=-i) 

ds' ds' 



(48) 



(49) 



which is a parity-odd but CP-even observable similar to the forward-backward asymmetry. The manifest expression 
for the longitudinal polarization asymmetry Pl in B q > — > Sll is derived as [87] 



Pl(s') 



2(1- ia) 1/2 
(1 + ^) as + nds 



Re 



C 



10" 



/+(<Z 2 



(50) 



It has been shown that this asymmetry is insensitive to the form factors in the masslcss limit for the lepton and can 
be approximated by 



Pl(s') 



2Re[C 9 e// C 10 ] 



|C 9 e// | 2 



0(C 7 ) 



1. 



(51) 



in view of the smallness of Wilson coefficient C-j compared with Cg and C\q . The distribution of the longitudinal 
polarization asymmetry Pj, in Bq — > Kq (1430)iZ as a function of q 2 are presented in Fig. 4, from which we indeed 
find that Pl(Bq — ► Kq (1430)e + e~) and Pl(Bq —> Kq (1430)/i + /i~) are close to —1 except the end points region. 

It is also useful to introduce the integrated longitudinal lepton polarization asymmetry (Apl) in order to characterize 
the typical value of longitudinal lepton polarization asymmetry 



{A PL ) = 



A PL {s')ds', 



(52) 
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FIG. 4: Lepton polarization asymmetries for Bo — » .R'o (1430)e + e , Bo — > .Ko (1430)/x + ^t and Bo — > Kg (1430)r + r as 
functions of squared momentum transfer q 2 based on light-cone QCD sum rules. 

with s' min = Amf /m 2 B and s' max = (m 2 B — m 2 s )/rn B . The numerical results of integrated longitudinal lepton polariza- 
tion asymmetry have been grouped in Table VIII, together with results in light-front quark model. From this table, 
we can observe that our results for (j4.pl) are in good agreement with that given by light-front quark model, which 
also indicates that this asymmetry is not sensitive to the decay form factors. 
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TABLE VIII: Numerical results of the integrated longitudinal lepton polarization asymmetry for Bo — * Kq ( 1430)ZT and B s — > 
/o(1500)i/ with I = e,n,T in the light-cone sum rules approach, , where the numbers estimated in light-front quark model [87] 
are also collected here. 



(Apl) 


Bo 


-> A' *(1430)e+e- 
-0.99 ±0.0 
-0.97 [87] 


B 3 


-> /o(1500)e + e" 
-0.99 ±0.0 




B 


-> /sT5(1430)At + M" 


B s - 


-+ /o(1500)/i + M~ 


(Apl) 




-0.96 ±0.0 




-0.96 ±0.0 






-0.95 [87] 






(Apl) 


Bo 


-> if *(1430)r+r- 
-0 03+ 00 
-0.03 [87] 


B s - 


-> /o(1500)r+r- 
-0.04 ±0.0 



VII. CONCLUSIONS 



Within the framework of light-cone sum rules, we analyze the form factors responsible for semi-lcptonic decays 
of B Q -> a (U50)Wi, B -> A" *(1430)/f, B s -> A^(1430)/P z and B s -> / (1500)Zf with I = e,fi,T up to the twist-3 
distribution amplitudes for the leading Fock state. Owing to the strong coupling of scalar mesons to the scalar current, 
the form factors associating with B — > S transition are approximately twice as large as that for the ones in the B — > P 
case. The form factors f+(q 2 ), f-{q 2 ) and /t(<7 2 ) calculated in this work verify the relations derived in the large recoil 
and heavy quark limit. 

Utilizing these form factors, we investigated the branching fractions of Bq — > ao(1450)/P;, Bq — ► Kq(1430)11, 
B s -> K*(U30)Wi and B s / (1500)/F. The magnitudes of BR(B Q -> a (1450)Z^) and -> A*(1430)Z^) 

can arrive of the order 10~ 4 , while the branching ratios of Bq — > Kq (1430)11 and B s — > /o(1500)ZZ are of the 
order 10~ 8 ~ 10~ 7 , which can marginally observed in the future experiments. The longitudinal lepton polarization 
asymmetries for Bq — > Kq(143Q)11 and B s — > /o(1500)ZZ are also considered in the SM. Our results for the asymmetry 
are in good agreement with that given by the light front quark model. The averaged asymmetries (Apl) for the 
final states including e + e~ and are almost equal to —1 except the end points region. However, the tau lepton 

polarization asymmetries are remarkably small and not measurable due to the efficiency for the detectability of the 
tauon. The theoretical predictions on the production properties presented in this work are very helpful to clarify the 
inner structures of scalar mesons as well as understanding the dynamics of strong interactions. 
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